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ABSTRACT 

We decipher intrinsic three-dimensional shape distributions of molecular clouds, cloud cores, 
Bok globules, and condensations using recently compiled catalogues of observed axis ratios for 
these objects mapped in carbon monoxide, ammonia, through optical selection, or in continuum 
dust emission. We apply statistical techniques to compare assumed intrinsic axis ratio distri¬ 
butions with observed projected axis ratio distributions. Intrinsically triaxial shapes produce 
projected distributions which agree with observations. Molecular clouds mapped in are 

intrinsically triaxial but more nearly prolate than oblate, while the smaller cloud cores, Bok 
globules, and condensations are also intrinsically triaxial but more nearly oblate than prolate. 


Subject headings: ISM: clouds — ISM: globules — ISM: structure —stars: formation 


1. Introduction 

Numerous catalogues have now compiled prop¬ 
erties of hundreds (or even thousands) of molec¬ 
ular clouds or cloud fragments with a large range 
of sizes, allowing meaningful statistical analysis of 
properties of molecular clouds, cloud cores, and 
smaller condensations. With such data sets, the 
distribution of apparent projected core axis ra¬ 
tio p can be used to constrain the intrinsic three- 
dimensional shapes. In a previous paper (Jones, 
Basu, & Dubinski 2001; hereafter Paper I), we 
investigated the shapes of molecular cloud cores 
mapped in NH 3 (Jijina, Myers, & Adams 1999) 
and cores mapped through optical selection (Lee 
& Myers 1999) via both statistical and analyti¬ 
cal methods. We found that strictly axisymmetric 
prolate or oblate shapes cores could not reproduce 
the observed projected axis ratios and that molec¬ 
ular cloud cores were triaxial. 

In this paper, we extend our previous work and 
conduct a statistical analysis of seven recent data 
sets which include a wide range of sizes of ob¬ 
jects, from molecular clouds with effective radius 


as large as 45 pc (Heyer, Carpenter, & Snell 2001) 
to submillimeter dust continuum maps of conden¬ 
sations with major axes as small as 2800 AU, or 
^ 0.01 pc (Motte et al. 2001). The information 
obtained about the intrinsic shapes of these ob¬ 
jects can yield insight into the physical processes 
which govern their evolution and subsequent star 
formation. We examine the catalogue of 

molecular regions in the outer Galaxy compiled 
by Heyer et al. (2001) in § 3.1, cores mapped in 
NH 3 and C^®0 in § 3.2 (Onishi et al. 1996; Tachi- 
hara, Mizuno, & Fukui 2000), catalogues of Bok 
globules (Clemens & Barvainis 1988; Bourke et al. 
1995) in § 3.3, and millimeter and submillimeter 
continuum maps of smaller protostellar condensa¬ 
tions (Motte, Andre, & Neri 1998; Motte et al. 
2001) in § 3.4. A discussion and summary are 
given in §§ 4 and 5. 
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2. Methods 

In general, a triaxial ellipsoid can be described 
by the equation 


where a is a constant and 1 > C ^ The 
geometrical analysis of Stark (1977) and Binney 
(1985) shows that such a body, when viewed in 
projection, has elliptical contours. Following Bin¬ 
ney (1985), the projection of a triaxial body when 
viewed from an observing angle {0, (p) (where the 
angles are defined on an imaginary viewing sphere 
and have their usual meaning in a spherical coor¬ 
dinate system) is found using the quantities 


. n , cos^d) 

siv?4) + 


B = cos6 sin2^ ^1 - 


U=y-^+cos~0Jj^. [i) 

The apparent axis ratio in projection then equals 


A + C-D 
A + C + D 


where D = ^{A — C)^ + B'^. Using these equa¬ 
tions, one can construct probability distributions 
for the projected axis ratio, assuming a large num¬ 
ber of randomly distributed viewing angles. 

We assign a Gaussian distribution of values for 
each axis ratio C and with a mean in the range 
[0,1], and standard deviation a typically equal to 
0 .1, consistent with our use of 10 bins to sample 
the data. We did test a range of cr from 0.05 to 
0 .2, and find that our conclusions do not change 
significantly within this range (see discussion in § 
4). The drawback to using a > 0.2 is that a rel¬ 
atively large fraction of the Gaussian distribution 
can fall outside the allowed range [0,1]. For exam¬ 
ple, a Gaussian distribution centered at 0.8 with a 
<j of 0.2 has « 16% of the data > 1 (see Paper I for 
an extended discussion of this issue). For similar 
reasons, we limit the cr = 0.1 analysis to the range 
of axis ratios [0.1, 0.9]. 


In order to find the best fit intrinsic distribu¬ 
tion of triaxial bodies, distributions of axis ratios 
with peak values and Co (for a given cr) are in¬ 
put into a Monte Carlo program. We typically 
employ at least lO'^ viewing angles to calculate 
the projected distribution for each individual pair 
of axis ratios, and at least lO'^ sets of axis ratios 
for each Gaussian distribution. This is more than 
sufficient for comparison with data sets sampled 
in 10 bins. The program produces the expected 
observed distribution which results from the as¬ 
sumed intrinsic distributions. We compare this 
output with the observed data sets via their 
values calculated by comparing the area of each 
bin to the area under the expected distribution at 
the location of each bin. Distributions of triaxial 
shapes for which the mean values Co = Co empha¬ 
size prolate objects, while those distributions with 
Co near 1 emphasize oblate objects. Furthermore, 
distributions with Co > ^(1 + Co) can be classi¬ 
fied as containing more nearly oblate than prolate 
objects. For greater details consult Paper I. 

3. Results 

3.1. Molecular Clouds 

Heyer et al. (2001) have catalogued the prop¬ 
erties of clouds and clumps using data from 

a survey of the outer Galaxy. The catalogue con¬ 
sists of 10 156 objects which include small, isolated 
clouds and clumps within larger clouds. Observa¬ 
tions from the outer Galaxy have the advantage 
that there is no distance ambiguity and the clouds 
are more widely spaced along the line of sight 
which eliminates problematic blending of emis¬ 
sion, allowing cloud properties to be determined 
more accurately. The FWHM beam size for this 
catalogue is 45 and there are an impressive 10 134 
objects with major and minor axes tabulated. 
Heyer et al. (2001) point out that the largest 
clouds, with effective radius re of about 10 pc 
or larger, are self-gravitating, whereas the smaller 
clouds (comprising the vast majority of observed 
objects) are not self-gravitating. This boundary 
also approximates the usual distinction between 
molecular clouds and giant molecular clouds, or 
GMC’s (see Blitz 1991; Williams, Blitz, & McKee 
2001). 

Given the evidence for a physical distinction 
based on size, we performed our analysis on the 
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entire data set as well as a subset with > 10 
pc, which corresponds to a mass range of approx¬ 
imately ~ 10^ to 10® Mq, typical of GMC’s. Fur¬ 
thermore, to avoid any pathological cases where 
an elliptical ht to the cloud shape can be a very 
poor approximation (M. H. Heyer, 2001, private 
communication), we restricted our sample to those 
clouds which span at least 10 spatial pixels in the 
observations (an even more restrictive threshold of 
20 pixels also yields essentially the same final re¬ 
sult). This criterion reduces the total set to 5685 
and the subset of GMC’s to 85. Our separate anal¬ 
yses can reveal whether there is any significant 
shape difference in the two populations. 

Figure 1 shows the results of the calcula¬ 
tions for the triaxial fitting of the GMC subset. 
The data set is best fit by distributions with axis 
ratios (^O)Co) = (0.2, 0.2) when a = 0.1. In order 
to determine whether values of (^o;Co) closer to 
zero would improve the Ht, we repeated the analy¬ 
sis with a value of cr = 0.05. (As the mean value of 
the Gaussian gets closer to the endpoints of the al¬ 
lowed range [0,1], a larger fraction of the Gaussian 
distribution falls outside this range, for a given cr. 
There is no ideal way to correct for this problem as 
explained in § 2 of this paper and in § 4 of Paper 
I). However, even with a = 0.05, the best fit mean 
axis ratios (^o, Co) are not closer to zero and agree 
with our result for ct = 0.1 within the estimated 
error. See § 4.1 for a discussion of the errors which 
we estimate to be a maximum of ± 0.1 in the mean 
value of each axis ratio. 

Figure 2 shows the result of the calculation 
for the complete set of axis ratios based on our 
selection criteria in the Heyer et al. (2001) cat¬ 
alogue. The best fit based on the x^ values is 
(Co,Co) = (0.3,0.3). 

For both the total set and for the subset based 
on large effective radius, the best-fit triaxial dis¬ 
tributions require = Co) as shown very clearly in 
Figure 1 and Figure 2. This means that the distri¬ 
butions emphasize prolate objects. However, since 
Co and Co are the means of distributions, most in¬ 
dividual objects cannot be considered strictly pro¬ 
late and are, in fact, triaxial. For the clouds with 
large re, the distributions which best fit the ob¬ 
servations have thinner objects (smaller C and C) 
than those that best fit the entire set. However, 
the difference is quite small and equal to our max¬ 
imum estimated error of ±0.1. 


Figure 3 compares the best fit distribution of 
p to the complete binned data set of Heyer et 
al. (2001). It also reveals that the histogram of 
projected shapes p of molecular clouds has some 
unique features. We recall that our previous anal¬ 
ysis (Paper I) of dense cores showed that an ob¬ 
served broad peak in the distribution at p ± 0.5 
and the presence of a significant number of objects 
near p = 1 favored triaxial, but more nearly oblate 
intrinsic shapes. Indeed, this pattern is reinforced 
in our subsequent study of other cores, Bok glob¬ 
ules, and protostellar condensations (see § 3.2 - 
3.4). However, the shapes of molecular clouds are 
distinct in that they have a very narrow peak, and 
at a low value p « 0.3. The narrow peak favors 
near-prolate objects, although a pure prolate cloud 
with ^ ^ = 0.3 yields a poor fit to the data, as 

shown in Figure 3. A pure prolate cloud would 
have too narrow a peak in the observed shape dis¬ 
tribution, as well as a higher probability than ob¬ 
served of a near-circular projection (see discussion 
and figures in § 2 of Paper I). In fact. Figure 3 
shows that even a Gaussian distribution of triax¬ 
ial objects implies a higher probability of detecting 
near-circular objects than observed. 

The cutoff in high values of p may be due to an 
actual cutoff in the distribution of intrinsic axis 
ratios (Co, Co) above some value, or could be due 
to some selection effect. One selection effect in the 
Heyer et al. (2001) sample is the fact that an ob¬ 
ject must span at least 5 pixels in the map to be 
classified as a cloud (additionally, we impose the 
higher threshold of 10 pixels for our shape anal¬ 
ysis). However, we see no evidence in the sam¬ 
ple for a trend toward greater circularity as clouds 
have smaller projected size. Another selection ef¬ 
fect is that the edges of the objects in the sample 
likely correspond to the the GO photodissociation 
boundary and not that of the H 2 gas (Heyer et al. 
2001). However, it is again not clear that this in 
any way biases against near-circular objects. 

3.2. Dense Cores 

Onishi et al. (1996) and Tachihara et al. (2000) 
have observed dense cloud cores in Taurus and 
Ophiuchus, respectively, in G®®0 with telescopes 
at Nagoya University. Since the same telescopes 
and technique were used to obtain these ratios, 
and both groups quote the axial ratios to a tenth 
of a parsec, we combined these two surveys in or- 
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Center of Gaussian 


Fig. 1.— Two-dimensional plot of inverse val- 
ues for triaxial core shape models applied to clouds 
with effective radius greater than 10 pc (and ob¬ 
served with at least 10 pixels) in the catalogue of 
Heyer et al. (2001). Note the symmetry of the 
figure about the line along which the two centers 
of the axis distributions are equal. For any point 
in the figure, the smaller axis ratio corresponds to 
^0 and the larger to Co- 
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Fig. 2.— Two-dimensional plot of inverse 
ues for triaxial core shape models applied to the 
entire catalogue of Heyer et al. (2001) based on 
our selection criterion. 


der to obtain a reasonable sized set for statistical 
analysis. The final set of data from these two re¬ 
gions includes 80 cores. The best fit determined 
by calculating the values is (Co, Co) = (0.4,0.9); 
see Figure 4. The best fit is compared with the 
binned data set in Figure 5. We performed this 
analysis with and without the sources for which 
there were embedded IR sources and obtained the 
same result. 

Although our best fit underestimates the ob¬ 
served axis ratios in the bin near p = 1 , it also 
overestimates the previous bin. The significant 
bin to bin variation is likely due to the relatively 
small size of this sample. However, the best fit 
agrees within the estimated errors, with the re¬ 
sults presented in Paper I using much larger sam¬ 
ples of NH 3 data (Jijina et al. 1999) and optical 
data (Lee & Myers 1999) for dense cores. Jijina 
et al. (1999) catalogued core properties for 264 
objects from NH3 observations and Lee & Myers 
(1999) catalogued properties for 406 dense cores 
from contour maps of optical extinction. The Ji¬ 
jina et al. data set is best fit by triaxial distri¬ 
butions with mean axis ratios (Co, Co) = (0.5,0.9), 
and the Lee & Myers data set is best fit by mean 
values (Co, Co) = (0.3,0.9). 

3.3. Bok Globules 

Bourke et al. (1995) catalogued physical char¬ 
acteristics of 169 isolated small molecular clouds 
(Bok globules) in NH3 from the southern sky. 
These observations were compiled to complement 
a similar study of 248 optically selected Bok glob¬ 
ules in the northern hemisphere by Clemens & 
Barvainis (1988). We note that Clemens & Bar- 
vainis (1988) also found that there was no correla¬ 
tion between the orientation of the clouds relative 
to the Galactic plane and their projected optical 
shapes. Ryden (1996) looked at the shapes of the 
globules from these two data sets based on the 
assumption of axisymmetry. She found that the 
Bok globules were consistent with oblate objects 
having an intrinsic mean axis ratio q = 0.3, or pro¬ 
late objects having q = 0.5. Ryden (1996) realized 
that these sets suffered from rounding errors. This 
affected the smallest major and minor axes, yield¬ 
ing an artificially large number of axis ratios near 
one. For example, Clemens & Barvainis (1988) 
measured the Palomar Observatory Sky Survey 
plates to compile their catalogue. They rounded 
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Fig. 3.— Comparison of the observed axis ratios 
from Heyer et al. (2001) with the best fit assum¬ 
ing triaxial clouds (solid line) and assuming a pure 
prolate cloud (dot-dashed line). The actual num¬ 
ber of observations in the bins is shown on the 
right-hand side vertical axis. 



Center of Gaussian 


Fig. 4.— Two-dimensional plot of inverse val¬ 
ues for triaxial core shape models applied to the 
combined surveys of Onishi et al. (1996) and 
Tachihara et al. (2000). 


the major and minor axes to the nearest millime¬ 
ter, which corresponds to an angular size of 1.12. 
The set is comprised of a relatively large number of 
small globules (42 out of 248) which have a major 
axis less than 2 mm. When these smallest values 
are rounded, the result is a large number of axis ra¬ 
tios erroneously equal to one or close to one. The 
top panel of Figure 6 shows very clearly that there 
are proportionally a large number number of ob¬ 
jects in the bin with the largest axis ratios. Obvi¬ 
ously, this problem does not affect the larger glob¬ 
ules as significantly. Ryden (1996) compensated 
for this problem by adding to each of the origi¬ 
nal major and minor axes a random error term A. 
She did this hundreds of times and obtained a new 
distribution of major and minor axes by taking an 
average from all her trials. See Ryden (1996) for 
additional details. We followed her procedure to 
correct for these rounding errors and obtained a 
new distribution of axis ratios for both data sets. 
Since the roundoff error is likely selected uniformly 
from within a fixed range, we used a random num¬ 
ber generator to select uniformly from a prescribed 
range of estimated error. For the Clemens & Bar- 
vainis (1988) and the Bourke et al. (1995) data 
sets, we used a range of error —0.56 < A < -|-0.56, 
and —0.25 < A < -1-0.25, respectively. Figure 6 
shows both the original binned data of axis ratios 
and the corrected distribution of axis ratios for the 
Clemens & Barvainis (1988) data. Notice that the 
large number of axis ratios near one in the origi¬ 
nal set has been dramatically reduced. We inves¬ 
tigated these two sets and found that intrinsic tri¬ 
axial shapes produced distributions which agreed 
with the observations. Figure 7 and Figure 8 show 
the results of the calculations with the triax¬ 
ial fitting for the Clemens & Barvainis (1988) and 
the Bourke et al. (1995) corrected data sets, re¬ 
spectively. The Clemens & Barvainis (1988) data 
set is best fit by distributions with mean axis ra¬ 
tios (^O)Co) = (0.4,0.9), and the Bourke et al. 
(1995) data set is also best fit by mean values 
(Co)Co) = (0.4, 0.9). The best fit for the Clemens 
& Barvainis (1988) data is also shown in Figure 6. 

3.4. Condensations Mapped in Millimeter- 
Submillimeter Wavelengths 

Motte et al. (1998) observed the p Ophiuchi 
main cloud at 1.3 mm with the IRAM 30 m 
telescope and were able to detect 58 individual 
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Fig. 5.— Comparison of the observed axis ratios 
from Onishi et al. (1996) and Tachihara et al. 
(2000) with the best fit assuming triaxial cores 
(solid line). The actual number of observations in 
the bins is shown on the right-hand side vertical 
axis. 
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Fig. 7.— Two-dimensional plot of inverse val- 
ues for triaxial core shape models applied to the 
catalogue of Clemens & Barvainis (1988) with the 
rounding errors corrected. 
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Fig. 6.— Clemens & Barvainis (1988) original 
binned normalized data set and the same data 
corrected for rounding errors. The original and 
corrected data are displayed in the top and bot¬ 
tom frames, respectively. The black solid line is 
the best fit from the triaxial distributions. The 
actual number of objects in the bins is shown on 
the right-hand side vertical axis. 
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Fig. 8.— Two-dimensional plot of inverse val¬ 
ues for triaxial core shape models applied to the 
catalogue of Bourke et al. (1995) with the round¬ 
ing errors corrected. 
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compact dusty objects with a fragmentation size 
scale of approximately 6000 AU. After we removed 
clumps which were thought to be composite in na¬ 
ture, there were 35 clumps for which projected ma¬ 
jor and minor axes were resolved. Although this 
is a small sample for statistical analysis, we still 
present the results of the calculations from the 
triaxial fitting in Figure 9. The Motte et al. (1998) 
data set is best fit by distributions with mean axis 
ratios (^o,Co) = (0.4,0.9). 

Motte et al. (2001) surveyed the protoclusters 
NGC 2068 and NGC 2071 in Orion B at 850 ^im 
and 450 pm with SCUBA (Submillimeter Com¬ 
mon User Bolometer Array) on the JCMT and 
were able to detect small objects (size ~ 5000 
AU) which they call condensations. There were 
64 condensations for which the projected major 
and minor axes were obtained. This is a remark¬ 
able data set which measures objects on size scales 
approaching that of our solar system. The re¬ 
sults of the calculations with the triaxial fit¬ 
ting are presented in Figure 10 for a = 0.1. The 
Motte et al. (2001) data set is best fit by distri¬ 
butions with mean axis ratios (CoiCo) = (0.4,0.9). 
Since the binned observed data of axis ratios does 
not have a strong peak, but instead looks rather 
flat, we repeated our analysis for a = 0.2. For 
cr = 0.2, the data set is best fit by distributions 
with mean axis ratios (^o,Co) = (0.5, 0.9). How¬ 
ever, the distribution with the mean value of 0.9 
and cr = 0.2, has 31% of the data values outside 
the allowed range. After the data is truncated, the 
remaining distribution will not have a Hnal mean 
at 0.9 and a = 0.2, and the remaining values in¬ 
side the allowed range must also be normalized. 
As a result, the normalized distribution is actu¬ 
ally peaked more than one would intuitively ex¬ 
pect. Consequently, the x^ values for a = 0.2 are 
no better than the values for cr = 0.1. The best 
fit for cr = 0.1 is compared with the observed data 
in Figure 11. 

Interestingly, Figure 11 shows that there are a 
number of objects with an apparent axis ratio near 
or exactly equal to one. A close examination of the 
original Motte et al. (2001) data reveals that the 
objects which fall into this bin are well above the 
resolution limit and that this large number is not 
a selection effect. Since this data set has a number 
of axis ratios near p = 1, it is interesting to specu¬ 
late whether or not intrinsic oblate objects (rather 


than near-oblate triaxial objects) may agree with 
the observations. We used two additional tests to 
investigate this hypothesis. First, we Ht the his¬ 
togram of observed shapes with a curve and carry 
out an analytic inversion to see if it is consistent 
with a distribution of intrinsic oblate shapes. The 
resulting intrinsic shape distribution tp^q) does go 
negative near <7 = 1 just as in the inversions per¬ 
formed on the dense core shape data of Jijina et al. 
(1999) and Lee & Myers (1999) in Paper I. This 
is due to the continuing presence of a decline in 
the observed axis ratio distribution (j){p) towards 
p = 1. To investigate further, we also experiment 
with an intrinsic axis distribution that is Gaus¬ 
sian in with ^0 = 0.4 and cr = 0.1, but has a 
fixed value ^ = 1 for all clouds, i.e., the clouds 
are all oblate with various degrees of flattening. 
In this case, our Monte Carlo program yields a 
probability distribution function which is plotted 
in Figure 11. The inverse x^ values for both the 
pure oblate objects and the best Ht triaxial ob¬ 
jects are « 0.2. The Motte et al. (2001) data, 
unlike the dense core data presented or reviewed 
in § 3.2, is close enough to being flat near p = 1 
and has a small enough sample of objects that we 
cannot entirely exclude the pure oblate hypothe¬ 
sis on the basis of the x^ test. If just one of the 
six objects (interestingly, one object has an em¬ 
bedded class 0 protostar) is removed from the bin 
closest to p = 1, the near-oblate triaxial objects 
provide a much better fit. Further observations 
and a larger sample will be necessary before the 
intrinsic shapes of these very dense objects can be 
determined with certainty. Nevertheless, we can 
conclude that molecular cloud condensations are 
either oblate or near-oblate in shape. 

4. Discussion 
4.1. Error Bounds 

In order to test the reliability of our results we 
recalculate the x^ values after randomly remov¬ 
ing 20% of the data and repeat this procedure 
ten times for each set. For every set investigated 
in this paper we obtain the same best fit mean 
axis ratios (■CojCo) for all the trials with the ex¬ 
ception of two trials with the Motte et al. (1998) 
data. This data set has only 35 non-composite 
objects for which the minor and major axis are 
published. When 20% of the data is randomly 
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Fig. 9.— Two-dimensional plot of inverse val¬ 
ues for triaxial core shape models applied to the 
survey of Motte et al. (1998) with the composite 
clumps removed. 



Fig. 10.— Two-dimensional plot of inverse val¬ 
ues for triaxial core shape models with cr = 0.1 
applied to the survey of Motte et al. (2001). 


removed, only 28 objects remain. Nevertheless, 
we obtain (^OjCo) = (0.4,0.9) for eight trials, 
(Co, Co) = (0.5,0.9) for one trial, and (Co, Co) = 
(0.4, 0.8) for the remaining trial. 

Additionally, adjusting the width of the Gaus¬ 
sian distributions in C and C within the range 
[0.05,0.2] yields at most a change of ±0.1 in the 
best fit values Co and Co- The largest variation is 
seen in the smallest data sets. The effect on the 
Motte et al. (2001) data is described in § 3.4. Us¬ 
ing cr = 0.05 allows us to test distributions with 
peaks near the boundaries 0 and 1, while a rel¬ 
atively large cr = 0.2 allows us to see if a better 
fit exists for the broadest observed distributions in 
projected axis ratio p. 

Altogether, our testing allows us to state an 
approximate maximum error in our best fit mean 
axis ratios of ±0.1. Interestingly, absolutely all of 
the best fits for the data sets of objects with size 
scale ~ 0.1 pc or smaller agree with one another 
within this range of error. This strengthens the 
case that these objects are triaxial but preferen¬ 
tially flattened in one direction and close to oblate¬ 
ness. A small degree of triaxiality seems necessary 
to explain the observed decline near p = 1, and the 
pure oblate hypothesis seems excluded for dense 
cores (Paper I and § 3.2). Due to the smaller 
number statistics and relatively large number of 
objects with p near one, the dense condensations 
are most likely to be compatible with pure oblate¬ 
ness. Future more extensive observations of this 
class of objects is necessary to settle the issue. 

Finally, there is the possibility that our results 
are biased by the fact that the spectral line or 
dust continuum data are probing emission from 
regions of varying opacity and/or temperature, so 
that the projected shape may not correspond ex¬ 
actly to the physical shape of the gas distribution. 
When data from other wavelengths become avail¬ 
able, it is possible our conclusions may change. 
This applies particularly to the molecular cloud 
data (§ 3.1), for which we have used a single cat¬ 
alogue. However, since we obtained very nearly 
the same result using a variety of tracers for the 
smaller cloud cores, Bok globules, and condensa¬ 
tions, the results in these cases seem robust. 
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4.2. Physical Implications 

This investigation shows that intrinsic triax- 
ial objects produce distributions which reasonably 
match observations of projected axis ratios for 
molecular clouds, molecular cloud cores, Bok glob¬ 
ules, and protostellar condensations. The results 
clearly fall into two categories: (1) on scales > 1 
pc, mapped in ^^CO, molecular clouds, includ¬ 
ing GMC’s, have triaxial shapes which are more 
closely prolate than oblate; (2) dense cores, Bok 
globules, and condensations, mapped in a vari¬ 
ety of tracers, on scales from few x 0.1 pc down 
to 0.01 pc, have triaxial shapes which are more 
closely oblate than prolate. The results about the 
latter objects reinforce our earlier finding (Paper 
I) from two other catalogues of dense core shapes. 
See Table 1 for a summary of the best fit axis 
ratios, {^o,Co), for each of the data sets we inves¬ 
tigated. The results from Paper I are included for 
comparison. 

The robust tendency for cores, Bok globules, 
and smaller condensations to have triaxial fits with 
^0 = 0.3 — 0.5, and = 0.9 implies that they 
are all preferentially flattened in one direction. 
This could be due to flattening along the direc¬ 
tion of a mean magnetic field, or due to signif¬ 
icant rotational support in the smallest objects. 
The magnetic field explanation for cores is attrac¬ 
tive as it implies that the observed near-alignment 
of core minor axes and magnetic held direction in 
Taurus (see Onishi et al. 1996) may be indica¬ 
tive of a more universal phenomenon. We also 
note that early submillimeter polarimetry of a few 
dense cores (Ward-Thompson et al. 2000) reveals 
a tendency toward alignment, but also a notice¬ 
able angular offset. This is interpreted as evidence 
for triaxiality of the cores (Basu 2000), which may 
still be preferentially battened along the direction 
of the magnetic held. 

While triaxiality is consistent with a nonequi¬ 
librium state, evolving due to external turbulence 
or internal gravity, the near-oblate shape also 
means that the objects may not be particularly 
far from equilibrium, and that oblate equilibrium 
models may act as a reasonable approximation to 
these objects. This can explain why the inter¬ 
nal structure of some Bok globules and pre-stellar 
cores can be closely or approximately ht by spher¬ 
ical equilibrium Bonnor-Ebert or near-equilibrium 
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oblate magnetic models (Alves, Lada, & Lada 
2001; Bacmann et al. 2000; Ciolek & Basu 2000; 
Zucconi, Walmsley, & Galli 2001). It is also con¬ 
sistent with the observed near-virial-equilibrium of 
most cores (Myers & Goodman 1988). 

We also note that the Bok globules, which are 
by dehnition isolated sites of star formation, have 
shapes that are not signihcantly diherent from 
that of molecular cloud cores and condensations 
embedded within larger clouds. This suggests that 
the environment in which the cores and condensa¬ 
tions are embedded plays a relatively insignihcant 
role in their dynamics, i.e., the external pressure 
from the parental cloud does not seem to be im¬ 
portant at this stage. 

The smallest objects in our study, the conden¬ 
sations mapped in millimeter and submillimeter 
continuum emission, may be the precursors to in¬ 
dividual stars since the mass spectrum appears 
to match the Initial Mass Function (IMF) com¬ 
piled by Salpeter (1955) over a certain mass range 
(Motte et al. 1998). The estimated triaxial but 
near-oblate shape of these objects are an impor¬ 
tant link in understanding the collapse process 
that leads to star formation. 

For the larger molecular cloud scale, we have 
utilized an exhaustive catalogue of the shapes of 
clouds in the outer Galaxy (Heyer et al. 2001). 
Although our study of molecular cloud shapes is 
based on this single available sample of projected 
axis ratios, and the result should be confirmed 
when other shape data become available, the sheer 
size of this catalogue is a strong point. The his¬ 
togram of observed axis ratios (Fig. 3) is very 
distinct from any of the other samples. It has a 
very sharp peak and a severe lack of objects with 
p > 0.5. While there may be some unknown selec¬ 
tion effect which biases against the observation of 
near-circular objects, we note that the the orien¬ 
tations of the projected shapes in the plane of the 
sky do appear to be truly random. Furthermore, 
we note that the earlier ^^GO catalogue of only 
23 clouds in Ophiuchus by Nozawa et al. (1991) 
that was utilized by Ryden (1996) has the same 
qualitative feature of a narrow peak near p = 0.3 
and a steep decline toward p = 1. 

Heyer et al. (2001) note that the vast majority 
of their clouds (all but the largest clouds, which 
we loosely label GMC’s) are not self-gravitating. 
They are either transient features or are held 
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Table 1 

Summary of Best Fit Mean Axis Ratios 


Data Set 

Object Type 

±0.1 

Co 

±0.1 

Heyer et al. (2001) 

clouds with effective radius > 10 pc 

0.2 

0.2 


complete set of clouds 

0.3 

0.3 

Onishi et al. (1996) & 

molecular cloud cores 

0.4 

0.9 

Tachihara et al. (2000) 

Jijina et al. (1999)®' 

molecular cloud cores 

0.5 

0.9 

Lee & Myers (1999)® 

molecular cloud cores 

0.3 

0.9 

Clemens & Barvainis (1998) 

Bok globules 

0.4 

0.9 

Bourke et al. (1995) 

Bok globules 

0.4 

0.9 

Motte et al. (1998) 

dense condensations 

0.4 

0.9 

Motte et al. (2001) 

dense condensations 

0.4 

0.9 


^Previous result from Paper I 



Fig. 11.— Comparison of the observed axis ra¬ 
tios from Motte et al. (2001) with the best fit 
assuming triaxial cores for cr = 0.1 (solid line) and 
assuming pure oblate cores (dashed line). The ac¬ 
tual number of observations in the bins is shown 
by the right-hand side vertical axis. 


together by external pressure. If these clouds 
are indeed brought together by large scale tur¬ 
bulence in the interstellar medium (or even con¬ 
fined for some time by an anisotropic ram pres¬ 
sure) we might expect that they have an elon¬ 
gated, filamentary shape (see e.g., Nagai, Inut- 
suka, & Miyama 1998; Balsara, Ward-Thompson, 
& Crutcher 2001; review by Shu et al. 1999). 
Since even the largest clouds seem to have these 
shapes, we surmise that all clouds may be brought 
together by external forcing (due to shock waves 
or turbulent motions for example), with only the 
largest clouds or densest regions within smaller 
clouds able to become self-gravitating. This ties 
in with the general picture of a rapid formation of 
molecular clouds due to external triggers (see e.g., 
Hartmann, Ballesteros-Paredes, & Bergin 2001; 
Pringle, Allen, & Lubow 2001). 

5. Summary 

Generally speaking, the observed decline in the 
observed probability distribution function (f>(j>) to¬ 
wards p = 1 favors triaxial rather than axisymmet- 
ric intrinsic shapes. In addition, objects observed 
on scales few x 0.1 pc and smaller have a broad 
peak and a significant number of objects observed 
with p « 1. This favors near-oblate triaxial ob- 
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jects, as shown in detail in Paper 1. In contrast, 
molecular clouds observed on scales > 1 pc have 
an observed (j){p) with a much sharper peak and 
a precipitous decline toward p = 1. This favors 
near-prolate triaxial objects. Reviews of the ex¬ 
pected distributions (j>{p) for various shapes can 
be found in Binney & Merrifield (1998) and Paper 
1. A summary of our best fits to the various data 
sets is given in Table 1. 

Our new results strengthen the finding that one 
of the best fit intrinsic axis ratios is always quite a 
bit larger than the other axis ratio for cores, con¬ 
densations, and Bok globules, which means that 
these objects are preferentially flattened in one di¬ 
rection, and close to oblateness. They may then 
not be particularly far removed from equilibrium 
or from oblate models often used to fit them. On 
the other hand, the much larger, lower density 
clouds have best fit distributions such that many 
objects are close to prolateness, consistent with 
the formation of these objects due to large scale 
external forcing. 

We thank John Dubinski for supplying the 
Monte Carlo program, as well as the anonymous 
referee and scientihc editor, Steven Shore whose 
comments and careful reading helped to improve 
the paper. This research was supported by a grant 
from NSERC, the Natural Sciences and Engineer¬ 
ing Research Council of Canada. C. E. J. acknowl¬ 
edges financial support from an NSERC postdoc¬ 
toral fellowship. 
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